Enhanced tracer transport by the spiral defect chaos state of a convecting fluid 
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To understand how spatiotemporal chaos may modify material transport, we use direct numerical 
simulations of the three-dimensional Boussinesq equations and of an advection-diftusion equation 
to study the transport of a passive tracer by the spiral defect chaos state of a convecting fluid. 
The simulations show that the transport is diffusive and is enhanced by the spatiotemporal chaos. 
The enhancement in tracer diffusivity follows two regimes. For large Peclet numbers (that is, small 
molecular diffusivities of the tracer), we find that the enhancement is proportional to the Peclet 
number. For small Peclet numbers, the enhancement is proportional to the square root of the 
Peclet number. We explain the presence of these two regimes in terms of how the local transport 
depends on the local wave numbers of the convection rolls. For large Peclet numbers, we further 
find that defects cause the tracer diffusivity to be enhanced locally in the direction orthogonal to 
the local wave vector but suppressed in the direction of the local wave vector. 

PACS numbers: 47.54.+r,47.27.Te,47.52.+j 



I. INTRODUCTION 



This paper addresses the transport of passive neutrally 
buoyant tracers in Rayleigh-Benard convection exhibit- 
ing spiral defect chaos — an example of spatiotemporal 
chaos that is characterized by disorder in both space and 
time Q 1^ 1^ . An important characteristic of such spa- 
tially disordered flows is that fluctuations in space play a 
significant role in their dynamics, resulting in advection 
of the passive tracers that is dependent in a complex fash- 
ion on space and time. The transport of passive tracers in 
such disordered flows is then governed by this advection 
in addition to molecular diffusion. The goal of this pa- 
per is to understand the net average transport of passive 
tracers as a function of the two competing mechanisms of 
advection by spatiotemporal chaos and molecular diffu- 
sion. Understanding material transport by spatiotempo- 
ral chaos is a problem that is of considerable importance 
in many branches of science and engineering. For exam- 
ple, an improved understanding may allow one to gain 
insight into heat and mass transport in atmospheric and 
oceanic flows and also in chemical engineering processes 
such as combustion. 

Previous studies of the properties of passive trans- 
port in convective flows have focused only on the steady 
and weakly oscillatory regimes. For example, in two- 
dimensional time-independent laminar Rayleigh-Benard 
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convection flow, experiments have shown that the trans- 
port is effectively diffusive in the long time limit, with 
an effective diffusivity that is greater than the molecular 
diffusivity by a factor that scales as the square root of 
the Peclet number (defined in Eq. (|14|l to be the ratio of 
the strength of advection to diffusion) (Jj . This enhance- 
ment, in the large Peclet number limit, has also been 
calculated theoretically by using the matched asymp- 
totic expansion method ^ .^J . In addition, higher-order 
corrections to the diffusion process, for arbitrary Peclet 
numbers, have been calculated numerically using the ho- 
mogenization method 0, • For nearly two-dimensional 
time-periodic convection, experiments near the onset of 
the oscillatory instability P| have shown that the trans- 
port is again effectively diffusive but with an effective 
diffusivity that depends linearly on the local amplitude 
of the roll oscillations ^3- This result has also been 
confirmed in theoretical work, which also identified the 
invariant structures of the flow that acted as templates 
for the motion of the tracers [ill IT^ . Passive tracer 
transport has also been studied in other types of laminar 
flow, including capillary waves generated by the Faraday 
instabihty 13,, 1 ^| and Taylor-Couette flow in a rotating 
annulus pL5l[la|. 

In this paper, the above transport studies are extended 
to flows that exhibit spatiotemporal chaos. It will be 
shown that the transport is globally diffusive and is en- 
hanced by the spatiotemporal chaos. However, unlike the 
case of laminar flows, the enhancement is found to follow 
two regimes. For large Peclet numbers (that is, small 
molecular diffusivities of the tracer), the enhancement 
is proportional to the Peclet number, whereas for small 



2 



Peclet numbers, the enhancement is proportional to the 
square root of the Peclet number. These two regimes are 
then explained by analyzing how the local transport de- 
pends on the local wave numbers of the convection rolls. 

The remainder of this paper is organized as follows. In 
Sect, m the equations governing Rayleigh-Benard con- 
vection and the transport of passive tracers are defined. 
In addition, direct numerical simulations of these equa- 
tions are discussed. In Sect. IIIII results from these sim- 
ulations are presented. In Sect. II V I conclusions are pre- 
sented. 



II. EQUATIONS AND ALGORITHMS 

A. Rayleigh-Benard Convection 

In a typical Rayleigh-Benard convection experiment, 
an incompressible fluid layer is confined between two hor- 
izontal plates, and is thermally driven far from equilib- 
rium by maintaining the bottom plate at a temperature 
that is higher than that of the top plate. As the tem- 
perature difference is increased, the fluid undergoes an 
instability to a state in which there is motion driven 
by the buoyancy forces. When the temperature differ- 
ence between the plates is above but near this convective 
threshold, a pattern comprising patches of locally parallel 
convection rolls forms with roll diameters that are close 
to the depth of the cell. When the temperature difference 
is increased, the fluid undergoes other instabilities that 
may result in the pattern developing an oscillatory or 
chaotic time-dependence. Finally when the temperature 
difference is increased further and if the aspect ratio F is 
larger than about 20 in boxes and 30 in cylinders, spiral 
defect chaos appears. This state is a disordered collec- 
tion of spirals that rotate in both directions and coexist 
with dynamical defects such as grain boundaries and dis- 
locations. Fig. n shows a numerically simulated instance 
of the spiral defect chaos state in a cylindrical geome- 
try. More generally, spiral defect chaos is an example of 
a kind of widely observed phenomenon called spatiotem- 
poral chaos that exhibits disorder in space and chaos in 
time. 

The evolution of the convecting fluid is governed to 
good approximation by the three-dimensional Boussinesq 
equations p^ . They are the combination of the incom- 
pressible Navier-Stokes and heat equations, with the fur- 
ther assumption that density variations are proportional 
to temperature variations and that this density variation 
appears only in the buoyancy force. Written in a dimen- 
sionless form, they are: 



a ^ [dt + u'V)u(a;,y, z,t) 
{dt + n'V)T{x,y,z,t) 
V'u 



0. 



RTziX) 
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FIG. 1: Example of spiral defect chaos observed in a numerical 
simulation described in Sect. Ill CI with insulating and no-slip 
boundaries, and with a spatial resolution of Aa; = 1/8 and a 
temporal resolution of At = 10"'^. The mid-plane tempera- 
ture field is plotted at time t — 500 for parameters e = 1.0, 
a — 1, and in a cylindrical geometry of aspect ratio F = 30. 
Dark regions correspond to cold sinking fluid, light regions to 
hot rising fluid. The spiral defect chaos planform is charac- 
terized by a disordered collection of spirals rotating in both 
directions and coexisting with dynamical defects such as grain 
boundaries and dislocations. 



fields, respectively. The variables x and y denote the 
horizontal coordinates, while the variable z denotes the 
vertical coordinate, with the unit vector z pointing in 
the direction opposite to the gravitational acceleration. 
The spatial units are measured in units of the cell depth 
d, and time is measured in units of the vertical thermal 
diffusion time <P / k, where k is the thermal diffusivity of 
the fiuid. The parameter R is the Rayleigh number, a 
dimensionless measure of the temperature difference AT 
across the top and bottom plates. 



R 



agd? 



AT, 



(4) 



where a is the thermal expansion coefficient, k the ther- 
mal diffusivity, and v the viscous diffusivity (kinematic 
viscosity) of the fluid. In this paper, the reduced Rayleigh 
number will also be frequently used. 



R ~ Rc 
Rc 



(5) 



The field u(a;, y, z, t) is the velocity field at point (x, y, z) 
at time t, while p and T are the pressure and temperature 



where Rc ~ 1708 is the critical Rayleigh number at the 
onset of convection in an infinite domain (ITj. The pa- 
rameter a is the Prandtl number, defined to be the ratio 
of the fluid's thermal to viscous diffusivities, 



(6) 
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The material walls are no-slip so that the velocity field 
satisfies 



u = 0, on all material walls. 



(7) 



The temperature field is constant on the top and bottom 
plates: 



T[x,y,z = T^,t\ =±i. 



(8) 



The lateral walls are assumed to be perfectly insulating, 
so that 



n • VT = 0, on lateral walls, 



(9) 



where n is the unit vector perpendicular to the lateral 
walls at a given point. The pressure field p has no asso- 
ciated boundary condition because it does not satisfy a 
dynamical equation. 

The influence of the lateral walls on the dynamics is 
determined by the dimensionless aspect ratio F, defined 
to be the half-width-to-depth ratio of the cell if it is rect- 
angular and the radius-to-depth ratio if it is cylindrical. 



B. Transport Equation 

The transport of passive neutrally buoyant tracers in a 
flow can be described by the advection-diffusion equation. 
Written in a dimensionless form, it is 



(10) 



The scalar field ip{x, y, z, t) is the passive tracer concen- 
tration at point (x, y, z) and time t. The velocity field u is 
obtained by solving the Boussinesq equations, Eqs. 

The parameter C is the Lewis number, defined to be 
the molecular diffusivity D of the tracers made dimen- 
sionless by the thermal diffusivity k of the fluid. 



D 

K 



(11) 



In this paper, small Lewis numbers in the range 10~^ < 
C < 10~^ will be considered. In comparison, the 
Lewis numbers of passive tracers used in previous con- 
vection experiments in water at approximately 300 K, 
namely, micrometer-sized latex spheres (vinyl toluene t- 
butylstyrene) and methylene blue dye, are C = 1.2x 10^^ 
and C = 3.9 x 10~^, respectively. 

The tracers are assumed to be passive, that is, their 
motions in the fluid do not modify the fluid's velocity 
field. The fiuid is also assumed to have negligible Soret 
and Dufour effects. The former refers to the additional 
passive tracer concentration current driven by gradients 
of the temperature field, whereas the latter refers to the 
additional heat current driven by gradients of the passive 
tracer concentration. In addition, the lateral walls are 
assumed to be impermeable to the tracers, so that 



where n is the unit vector perpendicular to the lateral 
walls at a given point. 

Eq. is also commonly written in the literature in 
an alternate but entirely equivalent form. By dividing 
it throughout by the product of a characteristic veloc- 
ity scale ||u|| and a characteristic length scale, Eq. H10|l 
becomes 



{di + u'V)-tJj{x,y,z,t) 



V 



(13) 



with t the rcscaled time, u the rescaled velocity field, and 
V the Peclet number, defined to be the dimensionless 
ratio of the relative importance of the advection of the 
tracers to their molecular diffusion. 



V 



c 



(14) 



(Note that the numerator in Eq. (|14|l contains a charac- 
teristic length scale — the depth of the cell — which is 
unity and is thus omitted.) 

Finally, it should also be noted that instead of study- 
ing the passive tracer concentration field -0 in the space 
coordinates defined in the laboratory frame (the Eule- 
rian approach), one could study the trajectories of each 
passive tracer individually (the Lagrangian approach) by 
integrating, for each passive tracer. 



dx(t) 
dt 



u[x(i),i]+77(t). 



(15) 



n • V'f/' = on lateral walls. 



(12) 



where x(t) is the position of the tracer [initially at x(0)], 
u is the Eulerian velocity field at space x(t) and time 
i, and ri{t) is a Langevin noise introduced to represent 
molecular diffusion. However, this approach is not pur- 
sued here because of the difficulties associated with inte- 
grating Eq. (|15|l . In fact, even if the velocity field u can 
be explicitly determined and has a very simple form, the 
tracer trajectories x can have very complicated dynamics 
0. 



C. Direct Numerical Simulations 

We use a parallel spectral element scheme to integrate 
the Boussinesq equations, Eqs. C[J~® , and the transport 
equation, Eq. H10|) . The scheme is second-order-accurate 
in time and is designed for rectangular, cylindrical, as 
well as more complex geometries with arbitrary lateral 
boundary conditions. Details of this scheme are available 
elsewhere 0| . For applications of this scheme to related 
problems in Rayleigh-Benard convection, see Refs. pol 

For small Lewis numbers C 1^ one well-known dif- 
ficulty associated with integrating Eq. ((Tn|l is 
that the spatial resolution Ax has to be very small. This 
scale is set by the smallest scale in the tracer concentra- 
tion field, such as the thickness of the interface where 
the tracer is initially zero on one side and unity on the 
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other. The interface is then stretched by a strain rate 
S ~ du/dx ~ ||u||, and its thickness is proportional to 
{ClSf^, that is, 



Aa 



1/2 

"uT 



C 



(16) 



For the spiral defect chaos states considered in this paper, 
the velocity magnitude, ||u||, is about 10. Thus, a simu- 
lation at the Lewis number of, say, C = 10"'^ will require 
Aa; ~ 10^^ in order to satisfy Eq. H16|) . However, cur- 
rent computational resources dictate that Aa; be about 
~ 10""'^, and in fact, for the results quoted in this paper, 
Aa; = 1/8. This problem is overcome by usinga filtering 
procedure developed by Fischer and Mullen "2711, and is 
described in Appendix^ Using this filter and maintain- 
ing Aa; — 1/8, a Lewis number as small ss, C — lO^'^ can 
be attained in a stable simulation. The accuracy of using 
the filter is also discussed in the Appendex. 



III. RESULTS 

The transport equation, Eq. (|1UI) . was integrated con- 
currently with the Boussinesq equations, Eqs. O-©, 
for the following parameters: the Rayleigh number var- 
ied from the onset of spiral defect chaos at i? « 3000 
to fully developed spiral defect chaos at i? w 4000, the 
Prandtl number a ~ \, and the Lewis number rang- 
ing from £ — 10^"^ to C — 10^^. The direct numeri- 
cal simulations have been performed in cylindrical three- 
dimensional cells of various aspect ratios. In this paper, 
data for an aspect ratio of F = 30 will be reported. The 
initial condition used for the passive tracer concentration 
field is a localized concentration at the center of the cell: 



'(/'(a;, y, z, < = 0) = exp 



y 



^2 1 



6A2 



(17) 



with A = 1/4 a small constant to ensure that the passive 
tracer concentration is initially localized. At time t = Q, 
the temperature, velocity, and pressure fields correspond 
to an asymptotic state of spiral defect chaos, that is, one 
that has been evolved from random thermal perturba- 
tions up to a time of OiV'^). In this paper, the focus will 
be on cells of large aspect ratio, F > 20. For these aspect 
ratios, the z-dependence of the passive tracer concentra- 
tion field was found to be essentially constant. As such, 
the z-dependence will be dropped in subsequent discus- 
sions and the passive tracer concentration ip(x,y,t) will 
be considered as a function of two-dimensional horizontal 
space and time. 

In Fig. 13 the evolution of the passive tracer concentra- 
tion field ip{x, y, t) at the mid-plane z = for the parame- 
ters R — 3500, <7 — \, and C — 10^^ is shown for various 
times t. The passive tracer concentration spreads out- 
ward with time in a non-uniform and non-axisymmetric 
way. In Sect. IIII Al this spreading is quantified globally 
by studying the mean square displacement of the passive 



R 
3074 
3500 
4270 



10"^ 

1.1 
1.0 
1.1 



lO""* 

1.1 
1.1 
1.1 



TABLE L The exponent 7 computed from the fit of the mean 
square displacement M2{t) to a power law ~ t'' for several 
different values of the Rayleigh number R and Lewis number 
£.. It is approximately unity in all instances. 



tracer concentration field. In Sect. 1111 lj| this spread- 
ing is shown to be characterized by normal diffusion. In 
Sect. IIII CI the local dependence of the spreading on the 
local wave number is discussed. 



A. Statistics of moments of passive tracer 
concentration 

The spreading of the passive tracers can be quantified 
by its mean square displacement M2{t), or the second 
moment, of the passive tracer concentration field, 



JO Jo 



{x){t)\^i;{r,0,t)rdrdd 



lo Iq^ ^(^' 



1 1) r dr dO 



(18) 



Here, x = (r, 9) is the polar coordinate with origin at the 
center of the cell. In practice, M2{t) is computed as the 
average over different instances (typically three to five) 
of V' obtained from different random initial conditions of 
spiral defect chaos, that is, the velocity, temperature, and 
pressure fields used ait — Q are different instances of fully 
developed spiral defect chaos. The quantity (x)(i) is the 
instantaneous center of mass of the tracer distribution, 



m = 



lo lo^ x'0(r, 9, t) r dr d9 
Jo Jo''Hr,9,t)rdrd9 ' 



(19) 



In Fig. 121 the mean square displacement M2 (t) is plot- 
ted for several different values of the Rayleigh number 
R and the Lewis numbers C — 1Q~^ and C — 10~^. It 
is found that, in all cases, the mean square displacement 
M2 {t) is directly proportional to the time t to a very good 
approximation. Least squares fits of M2{t) to power laws 
~ yield exponents 7 of approximately unity, as shown 
in Table n 

This implies that the spreading of the passive tracer 
concentration field can be described by a normal dif- 
fusive process. In other words, the averaged passive 
tracer concentration, ip{r,t), evolves according to the 
one-dimensional normal diffusion equation, 

dti^ir, t) = C*drri' (20) 

with C* an effective Lewis number that can be extracted 
from the mean square displacement as 

M2{t)=4:C*t. (21) 
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FIG. 2: Evolution of the passive tracer concentration field ^p{x,y, z — 0,t) for various times t, obtained by numerically solving 
Eq. mm . The Rayleigh number R — 3500, the Prandtl number a — 1, and the Lewis number C — 10~^, and the cylindrical 
cell has the aspect ratio F — 30. The initial condition for tp a.t t — is given by Eq. p7^ . 







10"^ 


10"^ 


10-1 


R 










3074 




0.29 


0.29 


0.38 


3500 




0.35 
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0.42 


0.41 


0.53 



TABLE IL The effective Lewis number C* computed from 
Eq. 12111 for various values of the Rayleigh number R and the 
Lewis number C. 



Several values of the effective Lewis number for various 
values of the parameters are tabulated in Table UTI 



The quantity 



A = 



C* -C 

c 



(22) 



is then a dimensionless measure of the enhancement in 
the molecular diffusivity of the passive tracer concentra- 
tion brought upon by the advection of the spiral defect 
chaotic flow. The goal of this paper can then be phrased 
as the calculation of how this enhancement varies as a 
function of both the properties of the advecting fluid (the 
Rayleigh number R and Prandtl number a) and the prop- 
erty of the passive tracer concentration [its Lewis num- 
ber, or equivalently, the Peclet number, P, c.f. Eq. (|14() ]. 

A^AiR,a,V). (23) 

The results of this calculation are shown in Fig. 01 
which depicts how the enhancement A varies vs. the 
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FIG. 3: The mean square displacement M2{t) of the passive 
tracer concentration field for two different Rayleigh numbers 
R = 3074 (top) and 3500 (bottom). The Prandtl number is 
CT = f in both cases. The triangle and circle symbols denote 
data for the Lewis number £, = 10~^ and £, = 10~^, respec- 
tively. The dashed lines have slopes of unity. The exponents 
obtained from power law fits of the data are given in Table Q] 
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FIG. 4: The dimensionless enhancement in molecular diffu- 
sivity A defined in Eq. vs. the Peclet number V for 

various Rayleigh numbers R and the Prandtl number a = \. 
Note that, when the Peclet number approaches zero (that is, 
when the advection becomes negligible so that the transport 
equation, Eq. IIIOII . is the diffusion equation), the enhance- 
ment should approach zero as well. Thus, the data points are 
expected to pass through the origin (P = 0, A = 0). 

Peclet number V for various Rayleigh numbers R and 
the Prandtl number a = 1. 

It is found that the enhancement A follows two dif- 
ferent scaling regimes in the Peclet number V. In the 
regime of large Peclet numbers, P~10^, the enhance- 
ment is found to scale linearly with the Peclet number, 

A cx P. (24) 



FIG. 5: The effective Lewis number £* — £ vs. the reduced 
Rayleigh number e. The solid line represents a power law with 
an exponent of 1/2. 

This implies that the effective Lewis number also scales 
linearly with the velocity magnitude of the flow, and is 
independent of the Lewis number £, 

£* oc ||u||. (25) 

[Eq. is easily obtained by dividing Eq. by the 
Lewis number £.] In addition, it is of interest to see 
how the effective Lewis number relates to the reduced 
Rayleigh number. This relation is plotted in Fig. [S] and 
it exhibits a square-root dependence, 

£*- roc 6^/2. (26) 

Thus, Eqs. I|25|) and H26I) together suggest that the char- 
acteristic velocity scale of spiral defect chaos scales with 
the reduced Rayleigh number as 

||u||ocei/2. (27) 

On the other hand, in the regime of small Peclet num- 
bers, P~10^, the enahancement is found to scale with 
the Peclet number as 

A cx (28) 

(Note that there are insufficient data to conclude whether 
the crossover between the two regimes is a continuous 
and gradual one or a discontinuous and sharp one.) 
This square root dependence of the enhancement on the 
Peclet number is similar to the result obtained exper- 
imentally ^ and calculated theoretically |^ in the 
spreading of passive tracers in time-independent convec- 
tion flows comprising straight parallel rolls. In this case, 
the enhancement can be attributed to the expulsion of 
the gradient of the passive tracer concentration from re- 
gions of closed stream lines 0. Near a separatrix be- 
tween two sets of closed stream lines, the only transport 
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of the passive tracers from one roll to the next comes 
from the random walks of the passive tracers that lie 
within a thin layer of width d of the roll boundary (the 
roll itself is of unit width in the dimensionless unit sys- 
tem adopted in this paper). Thus, a fraction d of pas- 
sive tracers contribute to an increase, dC, in the effective 
Lewis number of the diffusion. The width d can be es- 
timated from dimensional analysis 28] to be d^ ~ V~^. 
Combining these estimates leads immediately to Eq. (|28ll . 
Thus, the result above suggests that the gradient expul- 
sion mechanism near closed streamlines, although strictly 
derived in a time-independent convection flows compris- 
ing straight parallel rolls, may be a universal mechanism 
at sufficiently small Peclet numbers independent of the 
structure of the underlying flow field. 

In Sect, nil CI the origin of these two distinct regimes 
is discussed in terms of the dependence on the local 
wave number of the convection rolls. But before con- 
cluding this section, some details are presented in the 
way the least squares fits were performed. First, data 
from early times are ignored because of the presence 
of transients. One such transient effect could be that, 
at very early times prior to the turnover time scale 
Tc ^ ||u||~^ ~ 0{1Q~^), the passive tracers "feel" that 
they are being transported by a constant velocity field, 
and so will exhibit ballistic behavior with 7 = 2. There 
is then a crossover time in which 7 decreases to unity, 
and this regime is to be ignored too. Second, data from 
late times are also ignored because of finite size effects. 
The time at which finite size effects become important 
is chosen as the time at which the exponent 7, obtained 
from the logarithmic derivative 



lit) 



dl0g[M2it)] 

rflog(t) ' 



(29) 



deviates from approximately unity for a purely diffusive 
process with without advection (that is, whose diffusivity 
is chosen to match that of effective diffusivity of the above 
process). 



B. Normal diffusion vs. anomalous diffusion 

In this section, results are discussed from two other 
tests that show that the spreading process is indeed gov- 
erned by normal diffusion, and not anomalous diffusion. 
Diffusion is said to be anomalous when the mean square 
displacement is not proportional to time, that is, when 
M2{t) oc t"' with the exponent 7 1. Anomalous diffu- 
sion has been observed in the transport of passive trac- 
ers in cellular Taylor-Couette fiow in a rotating annu- 
lus 0, ! and in various other geophysical turbulent 
flows arising from the presence of Levy trajectories |29| |. 
However, results from this section show no evidence of 
anomalous diffusion in transport in spiral defect chaos 
for the range of Lewis numbers investigated. 

First, if the passive tracer concentration is spreading 
by normal diffusion and so obeys Eq. H2U|I , then it can be 
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FIG. 6: The scaled passive tracer concentration t^{r,t) vs. 
the scaled distance r'^ /At and scaled distance squared /At 
(inset). The various symbols denote different values of the 
times t. The data for times t = 30, t = 40, and t = IQ 
(triangles, squares, and circles, respectively) collapse onto the 
same curve, suggesting the validity of the Gaussian form in 
this time range. 



expressed in the form of a Gaussian, 



1 

- exp , 

t \AC*t 



(30) 



and a plot of the logarithm of the scaled passive tracer 
concentration log[tV'(r, i)] vs. the scaled distance i/r^/it 
for different times t will all collapse onto the same curve. 
This is indeed the shown in Fig. which shows 

the data at times t = 30, t = 40, and t = 70 (triangles, 
squares, and circles, respectively) collapsing onto the 
same straight line. However, data from an earlier time 
t — b (crosses) do not collapse onto the same straight 
line, presumably because of the presence of transient ef- 
fects. Similarly, data from a later time t = 100 (dots) do 
not collapse onto the same straight line, because of the 
presence of finite size effects. 

Second, possible deviations from the Gaussian behav- 
ior of the passive tracer concentration can be checked by 
looking at the higher-order moments. 



Jo Jo 1^ 



{y.){t)\'iil}{r,e,t)rdrde 



Jl^ C^j{r,e,t)rdrde 



(31) 



For normal diffusion, the higher-order moments scale like 

Mg{t) oc f^/^ , (32) 

and the ratio of this higher-order moment scaled to the 
second-order moment can be calculated to be 



Mg(t) 



2/q 



^l-3---{q-l) = {q-iy.\ 



(33) 



M2{t) 

which is a constant in time. In Fig.[7| this scaled ratio is 
plotted for g = 4, 6, and 8 as functions of time when the 
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FIG. 8: Illustrations showing the definition of the horizontal 
spreading orientation, O, at (a) Q = 7r/2 corresponding to 
spreading in the direction of the local wave vector k, and (b) 
= corresponding to spreading in the direction orthogonal 
to k. 



FIG. 7: The ratio of the scaled higher-order moments to the 
second-order moment, Mq{t)^^'' / Al2(t) , of the passive tracer 
concentration vs. time t for the parameters R — 3500, a = 1, 
and £ = 10~^, and for a purely diffusive process (dashed 
lines). From Eq. (PSJ, this ratio is 1.73 when g = 4, 2.47 
when q — 6, and 3.20 when q — 8. As q increases, the range 
of time for which the scaled moment stays constant in time 
decreases, because of finite size effects. 



Rayleigh number R = 3500, the Prandtl number cr = 1, 
and the Lewis number C = 10~^. The dashed hues show 
the corresponding quantity for a purely diffusive process 
with diffusivity chosen to match the former's effective dif- 
fusivity. (Because of finite size effects, the scaled higher- 
order moments unfortunately have only a small range for 
which they are constant in time. For example, at q = 8, 
this range is only 5~t~20.) The agreement of the two 
sets of data shows that, apart from finite size effects, 
there are no discernible deviations from Gaussian form 
for the passive tracer concentrations. 

Thus, both observations above suggest that the spread- 
ing of the passive tracer concentration is governed by 
normal diffusion. 



C. Wave number dependence of the passive tracer 
transport 

In this section, the existence of two different scaling 
regimes for the dimensionless enhancement in molecu- 
lar diffusivity, namely at large Peclet numbers given by 
Eq. H24() and at small Peclet numbers given by Eq. (|28|l . 
is investigated in terms of the local wave number depen- 
dence of the passive tracer concentration. To make this 
discussion more quantitative, first a quantity called the 
horizontal spreading orientation, Q{x,y), is calculated at 
every location in the cell. 



cos(e) 



(34) 



The subscript _L denotes the horizontal coordinates (x, y) 
and k(a;, y) is the local wave vector at location (x, y) in 
the planform |30j| . If the passive tracer concentration 
spreads in the direction of the local wave vector k, then, 
as illustrated in Fig. |S{a), the gradient V±ip will be or- 
thogonal to k, resulting in the local horizontal spreading 
orientation acquiring the value of = 7r/2. On the other 
hand, if the passive tracer concentration spreads in the 
direction orthogonal to k, then, as illustrated in Fig.|S{b), 
the local horizontal spreading orientation will be = 0. 
For a particular passive tracer concentration, the horizon- 
tal spreading orientation can be computed locally at ev- 
ery point in the mid-plane of the convection cell and then 
were sorted into bins to create a histogram. In Fig. |51 
such a distribution of the horizontal spreading orienta- 
tion, -P(0), is plotted for several values of the Lewis num- 
ber ranging from C — 10~^ to C = 10~^, the Rayleigh 
number R = 3500, and the Prandtl number cr = 1 at time 
t = 50. Consider first the distribution for the relatively 
large Lewis number of £ = lO"'^ (denoted by crosses). 
This distribution has a peak at = tt/2. This suggests 
that the spreading of the passive tracer concentration 
has the highest probability to be in the direction of the 
wave vector k [that is, as illustrated in the scenario of 
Fig. IHIa)] . This is consistent with the gradient expulsion 
mechanism near closed streamlines described earlier in 
Sect, nil Al However, at the smaller Lewis numbers of, 
say, C = 10~^, the distribution of the horizontal spread- 
ing orientation is visibly different. The distribution now 
has a peak at 8 = 0. In other words, the spreading of the 
passive tracer concentration has the highest probability 
in the direction orthogonal to the wave vector k [that 
is, as illustrated in scenario of Fig. (S^b)]. This shows 
that two different scaling regimes for the dimensionless 
enhancement in molecular diffusivity are associated with 
two different transport mechanisms. At large Lewis num- 
bers (or equivalently, small Peclet numbers), the trans- 
port is along the direction of the wave vector k by the 
gradient expulsion mechanism. At small Lewis numbers 
(or equivalently, large Peclet numbers), the transport is 
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FIG. 9: Distribution of horizontal spreading orientations 
P{0) for the Rayleigh number R = 3500, the Prandtl num- 
ber a — 1 and the Lewis number ranging from C = 10"^ to 

c = w-\ 



orthogonal to the wave vector k, presumably by advec- 
tion by the disordered flow. 

To determine why there is transport orthogonal to the 
wave vector k at small Lewis numbers (or cquivalently, 
large Peclet numbers), the following calculation was per- 
formed. The local wave numbers that correspond to loca- 
tions that exhibit spreading in the direction orthogonal to 
the wave vector k (that is, < < 77 where 77 = 0.01 is a 
small constant) were compared with those locations that 
exhibit spreading in the direction of the wave vector k 
(that is, 7r/2-77 < 6 < 7r/2). In Fig.HUl the distribution 
of wave numbers P{k, Q) for which the spreading occurs 
orthogonally to k (solid lines) is plotted together with 
the distribution for which spreading is along k (dashed 
lines) for a large Lewis number case [Fig. Iior al] and a 
small Lewis number case [Fig. llOf b)]. The key point to 
observe is that, in the small Lewis number case, there is 
a higher probability that spreading occurs orthogonal to 
the wave vector {& 0, solid line) than along the wave 
vector {Q 7r/2, dashed line) for wave numbers k « 1.5 
and k « 2.5. These wave numbers, far away from the 
mean wave number, correspond to to the occurrence of 
defects such as spiral cores, target cores, dislocations, 
etc. This suggests that the reason why gradient expul- 
sion ceases to be valid at large Peclet numbers is due to 
the enhanced transport of the passive tracers orthogonal 
to the local wave vectors by the defects in the pattern. 
However, the presence of defects is not sufficient to over- 
come the gradient expulsion mechanism at small Peclet 
numbers. 



IV. CONCLUSIONS 

In this paper, the spreading of a passive tracer concen- 
tration in a Rayleigh-Benard convection flow exhibiting 
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FIG. 10: Distribution of wave numbers P{k, O) for which 
spreading occurs orthogonally to the local wave vector k 
(solid lines) and in the direction of k (dashed lines), for (a) 
large Lewis number C = 10~^ and (b) small Lewis number 
C — 10~^ at the Rayleigh number R = 3500 and the Prandtl 
number a = 1. 



spiral defect chaos is studied for the first time. All previ- 
ous studies have dealt with time-independent or oscilla- 
tory flows. In the presence of advection by spiral defect 
chaos, we find that the spreading continues to be char- 
acterized by normal diffusion. The enhancement follows 
two regimes. When the Peclet number is large (that is, 
when the molecular diffusivity of the tracer is small), the 
enhancement is proportional to the Peclet number. This 
means that in the limit of large Peclet number the effec- 
tive diffusivity is independent of the molecular diffusiv- 
ity, and is proportional to the strength of the advection 
velocity field. When the Peclet numbers is small, the 
enhancement is proportional to the square root of the 
Peclet number. These results are explained in terms of 
the dependence of the transport on the local wave num- 
bers. It is found that tracers with small Peclet num- 
bers follow the gradient expulsion mechanism described 
previously in time-independent flows which predicts 
the square-root dependence. However, when the Peclet 
number becomes large, defects in the flow field became 
important and lead to enhanced transport orthogonal to 
the local wave vectors. 
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APPENDIX A: NUMERICAL DETAILS 

As mentioned in Sect. Ill ("I for small Lewis numbers 
£ ^ 1, a difficulty associated with integrating the trans- 
port equation, Eq. I|10|) , is that the spatial resolution Ax 
must be small. A larger Aa; can be used by employ- 
ing a simple filtering procedure developed by Fischer and 
Mullen [231 . At the end of each time step, a filter is ap- 
plied on an element-by-elenient basis to the passive tracer 
concentration field "ip- In one dimension, the filtered field 
can be written as 



,x 10 



F{^- a) = anN-iW + (1 - a)^, 



(Al) 



where the operator Hn-i first interpolates onto the 
mesh points for a polynomial of degree N —1 determined 
by the mesh spacing and then interpolates the result back 
onto the mesh points for a polynomial for degree N. In 
higher dimensions, the tensor product form of Eq. IjAip is 
used. Typical values of a used were in the range 0.05 < 
a < 0.2. This filtering procedure preserves inter-element 
continuity and spectral accuracy. Using this filter and 
maintaining Ax = 1/8, a stable Lewis number of up to 
£ = 10^^ could be attained. 

To verify that the filter allows the diffusion to be suffi- 
ciently and accurately resolved at small Lewis numbers, 
two sets of checks were performed. 

First, the passive tracer concentration field ip is in- 
spected for various values of a. In Fig. ^ a snapshot 
of the passive tracer concentration field ijj for C = 10^^ 
at a time i = 34 is shown for two values of a — 0.15 
(value used throughout in this paper) and a = 0.02. It 
can be seen that, although differences exist between the 
peak values of the filtered and the solution with a = 0.02, 
we can see that the extent and reach of the solution is 
largely unaffected it terms of the range of x it covers. 
Thus, when computing, say, the second moment of the 
passive tracer concentration field [Eq. l(TH|l ]. the results 
will be largely independent of a. This justifies that it is 
safe to use a filter parameter of as high as a = 0.15 when 
computing statistics of the passive scalar concentration 
field. 

Second, a local advection orientation, ^{x,y), is de- 
fined: 



cos(<I') 



Vip • u 



(A2) 



with u the velocity field. If the local passive tracer con- 
centration is being adverted by the local velocity and 
diffusion is not being sufficiently resolved, then the gra- 
dient of the former will be orthogonal to the local veloc- 
ity, and consequently, $ = 7r/2. On the other hand, if 
the local passive tracer concentration exhibits diffusion, 
then it will change in a direction perpendicular to the 




FIG. 11: Snapshot of the passive tracer concentration field 
ip{x,y = 0) for £ = 10~^ at time t — 34 for two different 
values of the filtering parameter a = 0.15 (solid line) and 
a = 0.02 (dashed line). 
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FIG. 12: Distribution of local advection orientations for var- 
ious Lewis numbers ranging from C = 10"* to C = 10"^ 



local velocity, yielding $ = 0. For small Lewis numbers 
where the effects of advection dominate over the effects 
of molecular diffusion, the distribution of the local ad- 
vection orientation, /'(<&), over the mid-plane of the cell, 
should exhibit a strong peak at <i> = 7r/2. This peak will 
then broaden as the Lewis number is increased, since the 
effects of diffusion cause the passive tracer concentration 
to spread out at all orientations relative to the local veloc- 
ity. The presence of this broadening in the distribution of 
the local advection orientation is then an indication that 
molecular diffusion has been sufficiently resolved. The 
distributions P(<&) for the various Lewis numbers rang- 
ing from C — 10"* to £ = 10~^ are plotted in Fig. ^1 
The distribution for C = 10~^ is distinctly different from 
that for C = 10~^, providing evidence that the molecular 
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diffusion at £ = 10""^ has been stably resolved, that is, 
that the chosen grid spacing A.x is sufficiently small for 
the simulation to be accurate. However, the relative sim- 
ilarity in the distributions for £ = 10~^ and £ = 10"'' 



suggests that diffusion for the latter case may not have 
been sufficiently resolved. Consequently, the smallest al- 
lowed Lewis number is set at £ = 10~^. 
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